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1.0  INTRODUCTION 


The  use  of  insensitive  munitions  (EM)  is  a  high  priority  among  the  military  services  from 
standpoints  of  both  safety  and  mission  effectiveness.  The  response  of  solid  rocket  motors  to 
bullet  and  fragment  impacts  is  an  important  factor  in  meeting  IM  propulsion  objectives. 
Considerable  effort  is  being  devoted  to  characterizing  the  bullet  and  fragment  vulnerability  of 
solid  rocket  motors,  and  to  developing  solid  rocket  motor  case  technologies  for  preventing  or 
lessening  the  violent  responses  of  rocket  motors  to  these  impact  loadings.  The  principle  hazards 
of  concern  include  prompt  shock-initiated  detonation,  delayed  detonation,  and  violent  explosion 
of  motor  propellants  from  bullet  and  fragment  impact.  Because  frill-scale  tests  are  costly,  fast¬ 
running  analytical  methods  are  required  to  characterize  the  response  of  solid  rocket  motors  to 
ballistic  impact  hazards  [1]. 

In  an  attempt  to  more  fully  understand  the  phenomena  involved  in  the  response  of  a  solid 
rocket  motor  to  bullet  and  fragment  impact,  a  first-principles-based  mathematical  model  was 
developed  to  determine  the  partitioning  of  the  kinetic  energy  of  the  impacting  projectile  among 
various  solid  rocket  motor  failure  modes.  This  kind  of  information  is  required  in  the  design  of  the 
various  components  of  IM  propulsion  systems,  such  as  the  motor  case  or  the  propellant  grain. 
Failure  modes  considered  in  the  analyses  included  case  perforation,  case  delamination,  and 
fragmentation  of  the  propellant  simulant  material.  Energies  involved  in  material  fragmentation 
were  calculated  using  a  fragmentation  scheme  based  on  a  procedure  developed  in  a  previous 
impact  study  utilizing  propellant  simulant  materials  [2-4],  These  quantities  were  computed  for  the 
fragmentations  caused  by  the  initial  projectile  impact  and  for  those  caused  by  the  impact  of  the 
residual  projectile  on  the  rocket  motor  rear  wall. 

Figure  1  shows  a  sketch  of  a  generic  rocket  motor  that  is  impacted  by  a  cylindrical 
projectile.  In  this  study,  we  consider  two  impact  scenarios.  In  the  first,  the  impact  velocity  is 
sufficiently  high  so  that  the  projectile  perforates  the  near  case  and  propellant  simulant  layers,  but 
not  high  enough  so  that  it  can  perforate  the  far  case  and  simulant  material  layers.  In  the  second. 
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the  impact  velocity  is  high  enough  so  that  the  projectile  is  able  to  perforate  both  sides  of  the 
rocket  motor  and  exit  with  some  residual  kinetic  energy.  In  this  manner,  the  model  developed 
herein  can  serve  as  the  basis  for  a  more  detailed  phenomenological  model  that  considers  a  much 
broader  range  of  impact  conditions.  However,  the  development  of  such  a  model  must  include  an 
experimental  component  whereby  empirical  evidence  is  obtained  to  guide  its  development. 
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2.0  SOLID  ROCKET  MOTOR  IMPACT  PHENOMENA 

Consider  the  impact  of  a  high  speed  projectile  on  a  generic  solid  rocket  motor.  The 
trajectory  of  the  projectile  is  such  that  the  projectile's  longitudinal  axis  is  perpendicular  to  that  of 
the  analog  motor  (i.e.  no  obliquity  and  no  yaw).  The  analog  motor  consists  of  a  layer  of 
propellant  simulant  material  bonded  to  a  filament  wound  case.  The  portion  of  the  filament  wound 
case  and  the  propellant  simulant  material  to  first  experience  the  effects  of  the  impact  are  referred 
to  as  the  'near*  case  and  simulant  material  layers;  corresponding  layers  to  feel  the  effects  of  the 
impact  subsequent  to  the  perforation  of  the  near  layers  are  referred  to  as  'far1  material  layers.  The 
following  discussion  of  the  phenomena  involved  in  the  perforation  and  fragmentation  of  solid 
rocket  motors  is  an  adaptation  of  the  discussion  presented  in  Ref.  [5]. 

Upon  impact,  strong  inelastic  waves  propagate  through  the  case  material.  These  waves 
propagate  radially  outward  from  the  impact  site  and  transversely  through  the  thickness  of  the 
solid  rocket  motor  wall.  The  radial  waves  damage  the  case  material  surrounding  the  impact  site 
This  damage  can  appear  as  delaminations,  fiber  pull-out,  fiber/matrix  fracture,  etc.  At  the 
case/simulant  interface,  a  portion  of  the  energy  associated  with  the  transverse  propagation  of  the 
waves  created  by  the  impact  is  reflected  back  into  the  case  material  in  the  form  of  release  waves 
while  the  remainder  of  the  energy  is  transmitted  into  the  simulant  material  in  the  form  of 
compression  waves.  These  waves  travel  through  the  propellant  simulant  material  and  eventually 
reflect  from  the  free  surface  as  tensile  waves.  The  backward  propagation  of  these  tensile  waves 
and  their  interaction  with  forward  propagating  compression  waves  causes  some  fracture  of  the 
rear  surface  of  the  simulant  material  to  occur. 

Simultaneous  to  the  motion  of  the  impact-induced  waves  in  the  case  and  simulant 
materials,  the  projectile  travels  through  the  case  and  simulant  materials  creating  a  tunnel  in  the 
simulant  material.  Its  velocity  naturally  decreases  as  some  of  its  energy  has  been  absorbed  by  the 
perforation  of  the  near  case  layer.  The  propellant  simulant  material  ahead  of  the  projectile  begins 
to  disintegrate  when  it  experiences  stresses  that  are  either  in  excess  of  its  ultimate  tensile  strength 
or  its  dynamic  fracture  toughness.  This  creates  a  cloud  of  propellant  simulant  debris  that  travels 


across  the  bore  and  impacts  the  far  layer  of  propellant  simulant  material.  If  sufficient  kinetic 
energy  is  available,  the  projectile  also  impacts  the  far  propellant  simulant  layer. 

The  impact  of  the  debris  cloud  on  the  far  propellant  simulant  material  layer  will  naturally 
damage  the  simulant  material.  If  the  kinetic  energy  of  the  original  projectile  is  sufficiently  high, 
then  it  is  entirely  likely  that  the  projectile  will  perforate  the  far  simulant  layer.  The  mode  of  far 
simulant  material  perforation  is  very  much  the  same  as  that  of  the  near  simulant  material  layer. 

The  projectile  tunnels  through  the  simulant  material,  eventually  disintegrating  the  plug  of  material 
that  forms  ahead  of  it.  Finally,  if  sufficient  kinetic  energy  remains,  the  projectile  may  perforate  the 
far  case  layer  and  exit  the  rocket  motor  with  some  residual  kinetic  energy. 

Based  on  this  discussion,  the  following  energy  balance  exists  for  the  initial  impact: 

Initial  Projectile  Kinetic  Energy  =  Energy  Required  to  Perforate  the  Near  Case  Layer 

+  Energy  Absorbed  by  Damage  to  Case  Material  Surrounding  the  Impact  Site 

+  Energy  Required  to  Fragment  the  Near  Propellant  Simulant  Material  Layer 

+  Energy  Required  for  Simulant  Debris  Motion 

+  Initial  Projectile  Residual  Kinetic  Energy  (if  any)  ( 1 ) 

If  there  is  non-zero  initial  residual  projectile  kinetic  energy,  then  it  is  absorbed  by  the  various 
phenomena  associated  with  the  impact  and  possible  perforation  of  the  far  material  layers.  Thus, 
the  following  energy  balance  can  be  written  for  the  final  residual  projectile  kinetic  energy: 

Final  Projectile  Residual  Kinetic  Energy  =  Initial  Projectile  Residual  Kinetic  Energy 

-  Energy  Required  to  Fragment  the  Far  Propellant  Simulant  Material  Layer 

-  Energy  Absorbed  by  Damage  to  Case  Material  Surrounding  the  Far  Impact  Site 

-  Energy  Required  to  Perforate  the  Far  Graphite/Epoxy  Case 

-  Energy  Required  for  Exterior  Simulant  Debris  Motion  After  Case  Perforation  (2) 

The  focus  of  the  remainder  of  this  paper  will  be  to  develop  mathematical  expressions  for  each  of 
the  terms  in  equations  (1)  and  (2).  However,  the  complexity  of  the  problem  dictates  that  certain 
assumptions  and  simplifications  are  required  to  bring  the  problem  down  to  solvable  proportions. 
These  assumptions  and  simplifications  are  discussed  in  the  next  section. 
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3.0  MODEL  ASSUMPTIONS  AND  SIMPLIFICATIONS 

In  developing  the  first-principles-based  model  of  kinetic  energy  partitioning  for  the  impact 
of  an  analog  solid  rocket  motor,  a  number  of  simplifying  assumptions  were  employed  These 
assumptions  are  discussed  in  the  following  paragraphs. 

First,  the  projectile  is  assumed  to  be  made  of  a  hardened  metallic  material  such  as  steel 
with  a  cylindrical  shape  and  a  blunt  nose.  The  filament  wound  case  is  presumed  to  consist  of 
graphite  fibers  wound  using  an  epoxy  as  the  matrix  material;  the  propellant  simulant  material  is 
presumed  to  be  essentially  rubber.  Further  comments  regarding  material  properties  will  be  made 
in  Chapter  6. 

Second,  the  cylindrical  analog  motor  is  modeled  as  two  parallel  flat  plates  (see  Figure  2) 
While  such  a  simplification  can  be  expected  to  have  some  influence  on  the  perforation  and 
fragmentation  response  characteristics,  the  deviations  of  the  flat  plate  characteristics  from  those  in 
which  a  cylindrical  geometry  were  used  are  expected  to  be  minimal  (see,  e  g.  [6]).  As  such,  it 
would  not  be  unreasonable  to  neglect  them  in  the  development  of  a  first-principles  based  model  of 
impact  response. 

Third,  the  projectile  is  assumed  to  remain  undeformed  as  it  perforates  the  near  filament 
wound  case  wall,  the  near  propellant  simulant  layer,  and,  in  the  event  that  it  exits  the  analog 
motor,  the  far  propellant  simulant  material  and  the  far  filament  wound  case  wall.  This  is  a 
reasonable  assumption  given  the  fact  that  the  impacting  fragment  is  presumed  to  be  made  of  a 
high-strength  hardened  metallic  material  such  as  steel,  whereas  the  filament  wound  case  and  the 
propellant  simulant  material  are  made  from  materials  of  comparable  or  lower  strength. 

Fourth,  the  hole  in  the  near,  and  if  appropriate  far,  case  wall  is  assumed  to  have  a  diameter 
equal  to  that  of  the  impacting  fragment.  This  is  also  a  reasonable  assumption  given  the  relative 
strengths  of  the  materials  involved.  In  addition,  if  the  impact  velocities  considered  sufficiently 
exceed  the  ballistic  limit  of  the  case  material,  but  not  to  the  point  where  the  materials  significantly 
melt  or  bum,  then  the  holes  created  in  the  case  walls  would  also  not  be  expected  to  significantly 
exceed  that  of  the  impacting  projectile. 


Fifth,  the  cavity  in  the  propellant  simulant  material  is  also  presumed  to  be  approximated  by 
the  shape  shown  in  Figure  3.  This  shape  is  based  on  photographs  obtained  from  experimental 
studies  in  which  blunt  nosed  projectile'-  impacted  slabs  of  propellant  simulant  material  The 
relationships  between  h,  hj,  and  6C  as  shown  in  Figure  3  will  be  discussed  in  the  next  section  As 
a  direct  corollary  of  the  second  assumption,  it  is  also  presumed  that  the  diameter  of  the  initial 
cylindrical  portion  of  the  cavity  shown  in  Figure  3  has  a  diameter  also  equal  to  that  of  the 
impacting  projectile. 

Sixth,  energy  absorbed  by  the  elastic  response  of  the  case  walls  following  the  impact  and 
that  lost  due  to  viscous  heat  dissipation  is  presumed  to  be  negligible  and  ignored.  The 
justification  for  the  neglect  of  energy  consumption  by  elastic  deformation  and  heat  dissipation  lies 
in  the  fact  that  the  impact  event  is  presumed  to  be  rapid  when  compared  with  the  strain  rates 
associated  with  elastic  response  and  with  heat  dissipation. 

Finally,  the  fragmented  far  simulant  material  is  presumed  to  remain  within  the  rocket 
motor,  and  in  the  event  of  far  case  wall  perforation,  only  the  projectile  is  ejected  from  the  motor 
This  appears  to  be  a  reasonable  assumption  because  the  far  graphite/epoxy  layer  confines  the 
simulant  material  and  generally  prevents  its  motion. 
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4.0  FRAGMENTATION  OF  PROPELLANT  SIMULANT  MATERIAL 

As  part  of  the  development  of  the  energy  partition  model,  a  model  of  the  fragmentation  of 
the  propellant  simulant  material  must  be  developed.  While  numerous  fragmentation  schemes 
currently  exist,  they  have  been  developed  primarily  to  model  the  processes  involved  in  the 
fragmentation  of  either  ductile  or  brittle  solid  materials  such  as  metals  and  rocks.  Thus,  the 
fragmentation  model  for  the  propellant  simulant  material  to  be  developed  here  is  semi-analytical, 
that  is,  it  is  based  on  empirical  evidence  but  presented  as  a  closed-form  analytical  expression. 

Based  on  the  results  of  a  curve-fitting  exercise  performed  using  the  data  provided  in  a 
recent  study  of  the  fragmentation  of  propellant  simulant  materials  [2-4],  we  postulate  that  the 
number  of  simulant  fragments  Nf(r)  with  a  radius  less  than  or  equal  to  a  specified  value  r  is  given 
as  follows: 

Nf(r)  =  A r-2  (3) 

where  A  is  a  constant  that  needs  to  be  determined  and  r  is  the  characteristic  radius  of  a  propellant 
fragment. 

To  determine  the  constant  A,  we  equate  the  total  mass  of  propellant  simulant  that  is 
fragmented  to  the  total  mass  of  fragments  whose  size  distribution  is  governed  by  equation  (3). 
Since  the  volume  of  propellant  simulant  fragments  is  equal  to  the  volume  of  the  cavity  created  in 
the  simulant  material,  and  since  the  cavity  may  be  said  to  consist  of  a  cylinder  resting  on  a  cone, 
we  have  the  following  equality  for  fragmented  simulant  mass: 

a  max  3  max 

7tR^(h-h1)/3  +  7tRph1  =  -  Jm(r)dNf (r)  =  -  J  m(r)-~dr  (4) 

3  min  3  min 

where  h  is  the  thickness  of  the  simulant  material  layer,  hj  is  the  depth  of  the  cylindrical  portion  of 
the  cavity,  Rp  is  the  radius  of  the  projectile,  and  Rc  is  the  radius  of  the  base  of  the  conical  portion 
of  the  cavity.  The  quantities  amax  and  amin  refer  to  the  radii  of  the  largest  and  smallest  simulant 
material  fragment;  m(r)=47tr^ppS/3  is  the  mass  of  a  fragment  with  radius  r  and  mass  density  ppS. 
Substituting  equation  (3)  into  equation  (4),  performing  the  integration,  and  solving  for  A  yields 
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(5) 


3\f  Rg(h-h,)/3  +  R^h 
8/1  amax-amin 

where  the  quantities  Rc,  hj,  amax,  and  any,,  are  still  to  be  determined. 

From  Figure  3,  is  simply  given  by 

Rc  *  (h-hj)tan0c  (6) 

where  0C  is  the  vertex  angle  of  the  conical  portion  of  the  simulant  material  cavity.  While  the 
value  of0c  is  undoubtedly  a  function  of  impact  velocity,  geometry,  and  material  properties,  a 
reasonable  approximation  for  0C  is  45  deg  [2-4]. 

The  functional  form  of  the  quantity  hj  is  also  based  on  the  empirical  evidence  provided  in 
Ref.  [2-4].  This  evidence  shows  that  hj~0  for  impact  velocities  greatly  in  excess  of  the  ballistic 
limit  of  the  simulant  material  layer  and  that  hjah  for  impact  velocities  near  the  ballistic  limit.  In 
other  words,  as  the  impact  velocity  is  increased  beyond  the  ballistic  limit  of  the  propellant 
simulant  material  layer,  the  depth  of  the  leading  cylindrical  portion  of  the  material  cavity  decreases 
from  a  value  that  is  slightly  less  than  the  full  thickness  of  the  layer  to  a  value  that  is  near  zero. 
While  there  are  any  number  of  functional  forms  that  meet  these  criteria,  the  form  chosen  for  use  in 
this  study  is  given  by 

h,/h  =  ea(l-V°/uBL>  (7) 

where  a  is  a  user-controlled  parameter  that  dictates  the  actual  shape  of  the  function,  VQ  is  an 
impact  velocity,  and  Ujjl  *s  the  ballistic  limit  velocity  of  the  simulant  material  layer.  The  actual 
definition  of  V0  and  the  means  of  calculating  the  ballistic  limit  velocity  of  the  propellant  simulant 
material  layer  UgL  sre  presented  in  the  next  section. 

While  there  is  clearly  no  theoretical  limit  for  the  size  of  a  smallest  fragment  created  during 
the  fragmentation  process  (after  all,  who  is  to  say  that  molecular-sized  particles  were  not 
created?),  common  sense  dictates  that  a  reasonable  cut-off  value  should  be  chosen.  For  the 
purposes  of  this  investigation,  a^n  was  arbitrarily  set  equal  to  0. 1  mm. 
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The  quantity  amax  can  be  obtained  by  assuming  that  there  is  only  one  particle  that  has  a 
radius  equal  to  amax.  Then,  by  the  definition  of  Nf,  we  have 

Nf(r  =  anMX)=l  (8) 

Substituting  for  Nf  according  to  equation  (3)  yields 


1  (9) 

a 

max 

so  that 

•max  “''A  (10) 

Substituting  for  A  according  to  equation  (S)  and  rearranging  yields 

aLx(amax-amin)-Tl  =  0  00 

where 


q  = 


(h-h,)/3+R2|l( 


(12) 


Equation  (1 1)  is  a  cubic  equation  for  amax.  To  solve  it,  we  first  rewrite  it  in  the  following  more 
standard  for: 


y3  +  py2  +  q  =  0 


(13) 


where  y^ma*,  p=-amjn,  and  q=-q.  Following  the  standard  procedure  for  the  solution  of  a  third 
degree  polynomial,  we  first  strive  to  eliminate  the  p  term.  This  is  accomplished  with  the 
substitution 


y  =  x-j  (14) 

in  equation  (13).  The  equation  obtained  is  then  written  as  follows: 

x3  +  3Hx  +  G  =  0  (15) 


<1  • 
mm 

9 


2p3  +  27q  2amin 

27  27 


(h-  h1)/3  +  Rph1 


(16a) 

(16b) 
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According  to  Cardan's  Formula  [7],  a  solution  to  equation  (1 5)  is  given  by 

1/3 

+ 

Since  equation  (IS)  has  three  roots,  it  is  also  desirous  to  determine  the  remaining  two  roots. 
However,  by  invoking  some  basic  principles  from  the  theory  of  equations,  we  can  infer  some 
information  regarding  the  nature  of  the  remaining  two  roots  without  having  to  solve  for  them. 
Consider  the  quantity 

D  =  -27(G2+4H3)  (18) 

IfD<0,  equation  (IS)  has  1  real  and  2  imaginary  roots;  if  D=0,  equation  (IS)  has  3  equal  real 
roots;  and,  if  DX),  then  equation  (IS)  has  3  distinct  real  roots.  Substituting  for  G  and  H 
according  to  equations  (16a,b)  into  equation  (18)  and  simplifying  results  in  the  following 
expression  for  D: 


(19) 


Since  and  q>0,  it  follows  that  D<0  and  equation  (15)  has  only  one  real  root.  Hence,  the 

other  two  roots  of  equation  (IS)  are  imaginary  and  can  be  ignored.  The  question  then  remains  as 
to  whether  the  root  given  by  equation  (15)  is  in  fact  the  real  root  or  one  the  imaginary  roots. 

Since  the  discriminant  in  each  of  the  terms  in  equation  (15)  is  simply  -D/27  and  since  D<0,  the 
discriminant  will  be  a  positive  number  so  that  a  non-imaginary  square  root  will  be  obtained. 

Hence,  the  root  given  in  equation  (15)  is  in  fact  the  sole  real  root  of  equation  ( 1 5).  The  value  of  x 
obtained  using  equation  (15)  can  then  be  used  to  define  amax  by  writing  equation  (14)  as  follows: 
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This  completes  the  development  of  the  fragmentation  scheme  used  in  this  study.  The  next 
section  addresses  the  partitioning  of  the  initial  kinetic  energy  of  the  impacting  projectile. 
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5.0  IMPACT  ENERGY  PARTITIONING 
In  this  section,  the  method  of  partitioning  the  kinetic  energy  of  the  impacting  projectiles  is 
discussed.  The  overall  approach  behind  the  method  is  to  base  the  analytical  formulation  on  the 
symbolic  equations  (1)  and  (2)  and,  in  doing  so,  incorporate  some  of  the  simplifying  assumptions 
discussed  previously.  Thus,  we  first  consider  equation  (1)  and  the  initial  impact  on  the  near  case 
and  simulant  material  layers,  and  follow  with  a  consideration  of  the  far  simulant  material  and  case 
layers. 


5.1  Impact  of  Near  Case  and  Simulant  Material  Laver 
In  this  part  of  the  development,  we  assume  that  the  residual  kinetic  energy  of  the  projectile 
after  passing  through  the  near  case  and  simulant  material  layers  is  not  sufficient  for  perforating  the 
far  simulant  material  and  case  layers.  Taking  into  account  the  simplifications  and  assumptions 
already  discussed,  equation  (1)  can  then  be  written  as 

(21) 


pinit  _  Fgr/ep  Fgr/ep  ,  pPs.cr  Fps,mv  Fresid 
proj  "  ^perf  +  ^delam  + frag  +  n  frag  +  ^proj 


where  EjJJJJj  is  the  initial  kinetic  energy  of  the  impacting  projectile  The  five  energies  on  the  right 

hand  side  of  equation  (21)  correspond  to  the  energies  required  to  perforate  the  near  case  layer,  to 
delaminate  the  case,  to  fragment  the  propellant  simulant  material,  to  move  the  debris  cloud 
consisting  of  simulant  material  fragments  across  the  bore,  and  the  residual  kinetic  energy  of  the 
projectile,  respectively. 

The  first  two  terms  on  the  right  hand  side  of  equation  (21)  are  simply  written  as 


^init  ~ 


~mpVn 

2  p  P 


<rfeP4mPVHL 


(22a) 


(22b) 


where  mp=7tRp2Lpppr  is  the  mass  of  the  impacting  projectile,  Vql is  the  ballistic  limit  of  the 
near  graphite/epoxy  case  layer,  and  ppr  is  the  mass  density  of  the  impacting  projectile. 
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To  determine  the  ballistic  limit  of  the  near  case  layer,  we  assume  that  the  propellant 
simulant  material  adjacent  to  the  near  graphite/epoxy  layer  does  not  significantly  affect  the  value 
of  the  ballistic  limit.  This  appears  to  be  a  reasonable  assumption  in  view  of  the  fact  that  the 
strength  of  the  propellant  simulant  material  is  significantly  less  than  that  of  the  case  material.  In 
deriving  an  expression  for  Vg^,  we  adopt  the  penetration  model  developed  by  Awerbach  [8] 
Although  this  model  was  originally  developed  for  metal  targets,  it  has  been  shown  to  yield 
acceptable  results  when  applied  to  composite  targets  as  well  [9].  In  this  model,  the  exit  velocity 
of  a  projectile  perforating  a  flat  target  plate  is  given  by 


Vf  = 


V.2  +  CT-»h 


V 


iI/2 


3p/2 


mo  /  pA 
^h  +  m0  /pA 


CTult 

3p/2 


(23) 


where  Vi  is  the  impact  velocity,  A  is  the  hole  diameter  in  the  target  plate,  mg  is  the  mass  of  the 
impacting  projectile,  and  h  is  the  thickness  of  the  target  plate;  p  and  oujt  are  the  density  of  the 
target  plate  and  ultimate  strength  of  the  target  plate  material.  In  the  event  that  V;  is  the  ballistic 
limit  velocity  of  the  target  plate,  then  Vf=0.  This  condition  can  be  used  to  solve  for  VgL  of  the 
graphite/epoxy  layer: 


VBL  = 


^  gult,gr/ep 
3pgr/ep  /  ^ 


mp  !  Pgr/epA 


\-3 


^gr/ep  +  mp  /  Pgr/epA  J 


1/2 


(24) 


where  hg,-/ep  is  the  thickness  of  the  graphite/epoxy  case.  Following  a  similar  line  of  reasoning, 
we  can  write  the  following  expression  for  the  ballistic  limit  of  the  propellant  simulant  material: 


CTult,ps  V  mp  /PpsA 
k  3ppS  /  2  J(^  h  +  mp  / ppsA  ^ 


(25) 


where  h  is  the  thickness  of  the  propellant  simulant  material.  Thus,  the  energy  required  to 
perforate  the  near  graphite/epoxy  layer  is  now  completely  defined. 
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The  energy  absorbed  by  the  delamination  of  the  case  material  is  a  bit  more  difficult  to 
define.  To  calculate  absorbed  energy,  we  must  first  be  able  to  calculate  the  actual  area  of 
delamination  that  occurs  in  the  case  material  as  a  function  of  velocity,  material  properties, 
geometry,  etc.  To  do  this  without  the  aid  of  empirical  data  is,  indeed,  quite  a  challenge. 

However,  not  to  be  deterred  from  our  ultimate  goal,  we  proceed  to  develop  a  means  for 
approximating  the  amount  of  delamination  that  occurs  that  should  be  suitable  for  a  first-principles 
based  analytical  model. 

Our  method  for  approximating  the  amount  of  delamination  is  based  on  a  trend  observed  in 
previous  experimental  studies  of  delamination  of  composite  materials  under  impact  loadings  and  a 
little  bit  of  common  sense.  It  has  been  previously  observed  that  as  impact  velocity  is  increased, 
the  total  delaminated  area  in  a  composite  material  increases  linearly  with  respect  to  the  kinetic 
energy  of  impact,  or  with  respect  to  the  square  of  the  impact  velocity  [10].  In  addition,  we  expect 
that  as  the  impact  velocity  greatly  exceeds  the  ballistic  limit  of  the  target  plate,  the  area  of 
delamination  will  approach  a  constant  value  (and  may  in  fact  begin  to  decrease  as  the  impact 
velocity  is  increased  still  more).  Thus,  the  total  area  of  delamination  can  be  approximated  as 

Adelam  =  Agfl^p)  (26) 

where  f[Vp)  is  a  function  with  the  anticipated  behavior;  A0  is  a  constant  denoting  the  asymptotic 
limit  of  the  function  defining  A^lam  An  elementary  function  that  exhibits  the  desired  behavior  is 

f(Vp)  =  l-e"PVp/VBL  (27) 

where  P  is  again  a  user-controlled  parameter  that  defines  the  shape  of  the  function.  Naturally,  a 
question  remains  as  to  the  value  of  Aq,  the  maximum  delaminated  area.  Experimental  studies 
have  shown  this  to  be  a  number  on  the  order  of  9  to  10  times  the  projected  area  of  the  impacting 
projectile  (see,  e  g.  [11]);  hence,  for  this  particular  investigation,  we  can  write 

Afci™  =  9Ap(l-e'PVVBL|  (28) 

However,  we  still  need  to  relate  A^lam t0  ^delam  ^is 's  accomplished  as  follows. 
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Suppose  the  near  and  far  case  layers  each  have  Nl  layers.  If  we  associate  a  critical 
surface  energy  density  yc  that  is  required  for  the  formation  of  a  delaminated  surface  and  assume 
that  the  area  of  delamination  is  constant  throughout  the  thickness  of  the  case,  then  we  can  write 
the  total  energy  required  for  the  delamination  of  the  near  case  layer  as 

Eg£--  2(N,.-l)Addamrc  (29) 

The  '2'  on  the  right  hand  side  of  equation  (29)  arises  from  the  fact  that  there  are  two  surfaces 
created  in  a  delamination  event.  Using  the  relationship  between  surface  energy  density  and  stress 
intensity  factor,  we  can  write 

Kric 

Yc=-nT  (30) 

where  c0  is  the  speed  of  sound  in  the  material  experiencing  fracture  (or,  as  in  this  case, 
delamination)  and  Knc  is  the  critical  mode  II  stress  intensity  factor.  Naturally,  there  are  various 
ways  of  defining  sound  speed  for  a  composite  material,  depending  on  the  direction  chosen,  etc. 
Since  this  is  a  first-principles  based  model,  it  will  be  assumed  that  the  sound  speed  can  be  written 
as 

c0  =  1 1  /  Pge  (31) 

where  Pge  is  the  density  of  the  graphite/epoxy  case  material  and  Ej  i  is  the  modulus  in  the  fiber 
direction.  Using  the  relationship  between  stress  intensity  factor  and  energy  release  rate,  the 
following  expression  is  written  for  the  mode  II  stress  intensity  factor  . 

KIIc  =  VEllGIIc  (32) 

where  Gflc  *s  the  critical  energy  release  rate  for  inter-layer  delamination. 

Thus,  substituting  equations  (28,30,3 1)  into  equation  (29)  results  in  the  following 
expression  for  the  amount  of  energy  absorbed  in  the  delamination  of  the  composite  case  material. 

9(NL-l)Ap(l-e-|iV'’/VBL]GIlc  (33) 
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To  derive  an  expression  for  E^r,  the  energy  required  to  fragment  the  propellant 

simulant  material,  we  follow  the  procedure  developed  in  Ref.  [3]  with  some  slight  modifications. 
Thus,  we  begin  by  writing 

amax  .vj 

EgJ>-  J  Es(r)^r  +  AcY  (34) 

3  min 

where  Es(r)  is  the  total  energy  consumption  associated  with  the  formation  of  a  single  propellant 
simulant  material  fragment  with  radius  r,  A<.  is  the  surface  area  of  the  simulant  material  cavity,  and 
Y  is  the  surface  energy  density  of  the  simulant  material.  In  Ref.  [3]  it  is  shown  that  Es  consists  of 
two  parts:  the  strain  energy  associated  with  the  deformation  of  a  potential  fragment  until  the  point 
of  rupture  and  the  surface  energy  associated  with  the  creation  of  the  fragment  surface.  Using  this 
formulation,  the  following  expression  for  Es(r)  is  obtained: 

Es(r)  =  Y*YT2  (35) 

where  y,  as  used  in  equations  (34)  and  (35),  is  written  as 

y  =  — —  (36) 

2ppscs,ps 

where  Kn>  pps,  and  cSjpS=V(Kps/pps)  are  the  mode  I  dynamic  fracture  toughness,  density,  and 
sound  speed  of  the  propellant  simulant  material;  Kps  is  the  bulk  modulus  of  the  propellant 
simulant  material.  Substituting  for  Nf(r)  according  to  equations  (3)  and  (5)  and  performing  the 

required  manipulations  results  in  the  following  expression  for  Ej^: 


Pps,cr  _ 
frag  “ 


K 


ro.ps 


2pPScs,ps  J 


5x 


Rc  (h  ~  ht)/  3  +  Rp 

amax  ~  amin 


taf  lOE-V  *Rc  jRj+ai-h,)1 

V  a  min  J 


(37) 
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Since  the  fragmentation  of  the  propellant  simulant  material  is  due  to  the  impact  of  the 
projectile  after  it  perforates  the  graphite/epoxy  case,  the  'impact  velocity*  used  in  equation  (37)  to 
calculate  h|  is  Vf,  the  reduced  velocity  oi  the  projectile  as  defined  by  equation  (23). 

This  completes  the  derivation  of  the  expression  for  the  energy  absorbed  in  the 
fragmentation  of  the  propellant  simulant  material.  The  final  expression  to  be  derived  is  that  for 
the  kinetic  energy  of  the  moving  debris  cloud.  This  energy  expression  is  written  as  follows: 


pps.mv 
*  fiag 


(38) 


where  m(|c  is  the  total  mass  of  the  propellant  simulant  fragments  in  the  debris  cloud  and  V(jc  is 
the  velocity  of  the  debris  cloud  center  of  mass.  Based  on  the  development  thus  far,  the  mass  of 
the  debris  cloud  material  can  be  written  simply  as 

mdc  =Pps7tRc^h~hi)/3  +  pPsRphl  (39) 


The  velocity  of  the  debris  cloud  is  also  obtained  by  following  the  procedure  in  Ref.  [3],  First,  we 
consider  the  following  energy  balance: 


(40) 


where  Eavaji  is  the  amount  of  energy  delivered  to  the  fragmentation  volume,  i.e.  it  is  the  amount 
of  energy  available  for  fragmenting  the  propellant  simulant  material  and  its  motion: 


E 


avail 


=  Eps’cr 
11  frag 


+  E 


ps,mv 

frag 


(41) 


Thus,  if  we  can  derive  an  independent  expression  for  Eavaji,  we  can  directly  solve  for  Vjc  by 
rewriting  equation  (41)  as  follows: 


mdc^dc  -  Eavail  Efra 


,cr 
frag 


(42) 


where  E£f 


is  given  by  equation  (37).  Following  the  procedure  in  Ref. 


[3],  we  write 


Eavai!  4mS<2Vr’>: 


(43) 
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where  m^  is  the  mass  of  the  compressed  part  of  the  propellant  simulant  material  and  V™11  is  the 
ps  ps 

velocity  imparted  to  the  propellant  simulant  material,  that  is,  it  is  the  particle  velocity  in  the 

simulant  material  induced  by  the  passage  of  the  initial  shock  wave  through  the  material.  The 

presence  of  the  '2'  in  front  of  the  velocity  term  in  equation  (43)  is  due  to  the  fact  that  as  the  shock 

wave  in  the  propellant  simulant  material  travels  through  the  material  and  reflects  from  the  rear 

free  surface,  it  induces  motion  in  the  simulant  material  that  is  characterized  by  a  velocity  equal  to 

twice  the  particle  velocity  in  the  material. 

Using  the  three  shock  jump  conditions  to  characterize  the  impact  of  the  projectile  on  the 
propellant  simulant  material  after  perforation  of  the  graphite/epoxy  case  and  continuity  of 
pressure  and  material  at  the  impact  interface,  we  obtain  the  following  expression  for  particle 
velocity  in  the  propellant  simulant  material: 


where 


(45a) 


b  =  2kprVf  +cSipr  +cS(PS(^i)  (45b) 

Ppr 

c  =  cs,pr^f  +  kpr^f  (45c) 

where  the  k  constants  are  the  slopes  of  the  linear  shock  wave  speed-particle  velocity  re.ationships 
for  the  projectile  and  propellant  simulant  materials  and  the  cs  =  ^K/p  are  the  adiabatic  sound 

speeds  of  the  materials. 

The  mass  of  the  compressed  portion  of  the  propellant  simulant  material  can  be 
approximated  as  the  mass  of  the  material  directly  beneath  the  projectile  as  it  penetrates  into  the 
simulant  material,  that  is. 
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n,s  =  e7tRphpps 


where  e  is  another  user-controlled  parameter.  Substituting  equations  (43,44,46)  into  equation 
(42)  yields  the  following  expression  for  Vfc. 


V*.  = 


f^PpsK^-E*;, 


mdc  1 2 


Finally,  the  expression  for  the  residual  velocity  of  the  projectile  is  written  as 


= 


mp  /  2 


where  E^d  is  the  energy  remaining  after  the  initial  impact  energy  has  been  partitioned  among 

the  energy  absorbing  processes  discussed,  that  is, 

presid  _pinit  /pgr/ep  ,  pgr/ep  ,  pPS.cr  pps.rnv, 

Eproj  -Eproj  (fcperf  +  Edelam +  Efrag  +bfxag  > 


5.2  Impact  of  Far  Case  and  Simulant  Material  Laver 


If  the  impact  velocity  is  significantly  larger  than  the  ballistic  limit  of  the  graphite/epoxy 
case,  then  the  possibility  exists  that  the  projectile  will  exit  the  rear  of  the  solid  rocket  motor.  In 
this  case,  equation  (2)  can  be  rewritten  as 

Eresid  =  £ps,cr  Egr/ep  ggr/ep  gfinal  (50) 

proj  frag  deiam  perf  proj  v  ’ 

In  writing  equation  (50),  it  is  noted  that  there  is  no  energy  allocated  to  the  motion  of  the 

propellant  simulant  fragments  created  by  the  far  impact.  This  is  due  to  the  presence  of  the  far 

case  layer  which  prevents  the  motion  of  the  far  simulant  material  fragments.  The  quantities 

Eres’d  E817!*5,  and  E?*’cr  are  defined  by  equations  (49),  (22b),  and  (37),  respectively.  However, 
proj  pen  *  rag 
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it  is  noted  that  in  evaluating  e££“  for  the  far  layer  using  equation  (37),  the  quantity  h\  is 

evaluated  using  an  'impact1  velocity  of  Vres,  which  is  the  residual  velocity  of  the  projectile  after 
exiting  the  near  case  and  propellant  simulant  material  layers. 

The  next  energy  absorbing  mechanism  considered  is  the  delamination  of  the  far 
graphite/epoxy  case  material.  Following  the  procedure  used  previously,  we  arrive  at  the 
expression 

E&  =  9<NL  -1)Ap(l-e_pV,/VBL  )Gnc  (51) 

where  V  has  been  written  instead  of  Vp  to  take  into  account  the  reduced  velocity  of  the 
projectile.  This  reduction  in  velocity  is  due  to  the  energy  loss  in  fragmenting  the  far  simulant 
material  layer  and  is  obtained  from  the  following  energy  balance: 


|rap(v') 
Substituting  for  E 


2  _  presid  _  pl»,cr 
proj  Cfrag 

j^d  according  to  equation  (48)  and  solving  for  V  yields 


(52) 


pps.cr 

irag 

nip  /  2 


(53) 


where  E^r  is  given  by  equation  (37)  with  hj  again  being  evaluated  using  an  'impact'  velocity  of 

Vres.  Following  the  calculation  of  V,  it  is  worthwhile  to  check  and  see  if  V>Vbl  If  not,  there 
is  not  enough  energy  remaining  to  perforate  the  far  case  layer  and  the  projectile  becomes 
embedded  in  the  far  propellant  simulant  layer.  If  in  fact  it  is  true  that  V^Vbl-  lhen  it  is  indeed 

a 

likely  that  the  projectile  will  perforate  the  far  case  layer  and  exit  the  solid  rocket  motor. 

The  final  kinetic  energy  of  the  projectile  as  it  exits  the  solid  rocket  motor  and  its  final  exit 
velocity  are  then  readily  obtained  using 


and 
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respectively.  This  completes  the  analytical  development  of  the  energy  partitioning  model.  The 
next  section  presents  some  numerical  results  using  the  equations  derived  herein. 


6.0  NUMERICAL  RESULTS  AND  DISCUSSION 

A  FORTRAN  program  was  written  to  implement  the  equations  developed  in  the  previous 
section;  its  listing  is  provided  in  Appendix  A.  Typical  output  files  are  provided  in  Appendix  B.  In 
the  implementation  of  the  energy  partition  model  developed  herein,  the  following  values  were 
used  for  the  user-controlled  constants  a,  (3,  and  e. 

a  =  0.25,  (3=1.0,  e  =  0.5  (56) 

Naturally,  experimental  information  is  required  to  ultimately  justify  these  choices  or  dictate  the 
use  of  other  values.  These  specific  values  were  chosen  because  their  use  allowed  the  predictions 
of  the  model  to  be  consistent  with  expected  trends  in  the  output  as  well  as  with  the  laws  of 
physics.  For  example,  the  specified  choice  of  e  ensured  that  the  velocity  of  the  projectile 
decreased  monotonically  as  the  projectile  moves  through  the  rocket  motor  An  improper  choice 
of  e  could  actually  result  in  an  increase  in  the  projectile  velocity  after  perforating  the  near  case 
layer  and  fragmenting  the  near  propellant  simulant  material.  In  addition,  these  choices  resulted  in 
ballistic  limit  values  for  the  graphite/epoxy  and  the  propellant  simulant  that  are  consistent  with 
those  found  in  the  literature.  Specifically,  the  model  developed  herein  predicted  that  Vql=250 
m/sec  and  Uql»50  m/sec.  These  values  are  consistent  with  those  cited  in  Refs  [9]  and  [2-4], 
respectively. 

The  material  properties  used  in  the  numerical  simulations  are  given  in  Table  1 .  The 
material  properties  were  obtained  by  considering  the  fiber  material  to  be  EM6  graphite  while  the 
matrix  material  was  taken  to  be  3501-6  epoxy.  The  material  property  values  presented  in  this 
Table  are  representative  values  obtained  by  averaging  several  values  obtained  during  a  literature 
survey  (see,  e  g.,  [2-4,12-27]).  The  properties  in  the  column  next  to  the  propellant  simulant 
material  represent  typical  thermodynamic  and  mechanical  properties  of  actual  propellant  materials. 
Worthy  of  note  in  Table  1  is  the  fact  that,  for  the  most  part,  the  material  property  values  of  the 
propellant  simulant  material  closely  resemble  those  of  an  actual  generic  propellant  material. 
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Table 

.  Material  Properties 

Property 

Units 

Propellant 

Simulant 

1  ■ 

isjsm 

km/sec 

3.58 

2.21 

2.18 

k 

0.91 

1.87 

2.11 

rft 

0.50 

1.44 

1.15 

P 

1.53 

1.78 

1.73 

E 

GPa 

177.5 

2.40 

498 

v 

,  ,r 

0.39 

0.45 

0.40 

K 

GPa 

N/A 

8.17 

8.31 

G 

GPa 

4.62 

0.81 

1.78 

mm 

MPaVm 

N/A 

N/A 

1.28 

J/m2 

400 

N/A 

N/A 

_ _ 

MPa 

N/A 

0.37 

N/A 

mm 

MPa 

1.38 

0.59 

N/A 

Figure  4  shows  plots  of  the  propellant  simulant  debris  cloud  velocity  and  the  final 
projectile  velocity  as  functions  of  the  initial  impact  velocity.  The  relationship  between  the  debris 
cloud  velocity  and  the  impact  velocity  is  seen  to  be  relatively  constant  until  the  impact  speed  at 
which  the  projectile  is  able  to  perforate  the  far  case  layer  as  well  (approximately  420  m/sec)  and 
then  becomes  nearly  linear.  The  plot  of  the  final  projectile  velocity  is  also  seen  to  change  near  an 
impact  speed  of  420  m/sec.  However,  the  change  in  the  plot  of  the  final  projectile  velocity  is 
more  dramatic:  there  is  a  pronounced  'dip'  in  the  curve.  This  odd  behavior  can  be  explained  by 
the  following  considerations. 

Until  the  projectile  perforates  the  far  case  layer,  its  residual  velocity  steadily  increases 
from  its  near  zero  value  when  the  projectile  just  barely  perforates  the  near  case  layer.  In  this 
initial  regime,  the  projectile's  residual  velocity  is  below  the  ballistic  limit  of  the  case  material.  As 
such,  the  projectile  remains  embedded  in  the  far  simulant  layer  without  perforating  the  far  case 
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layer.  Naturally,  as  the  impact  velocity  is  increased,  the  depth  of  embeddment  also  increases. 
Finally,  when  the  impact  velocity  is  high  enough  so  that  the  residual  velocity  exceeds  the  far  case 
material's  ballistic  limit,  the  projectile  perforates  the  far  case  layer.  However,  in  doing  so  a 
substantial  amount  of  energy  is  absorbed  by  the  perforation  process.  Hence,  at  impact  velocities 
just  above  those  required  for  far  case  layer  perforation  the  final  projectile  velocity  will  in  fact  be 
relatively  low.  As  the  impact  velocity  continues  to  increase  so  does  the  final  projectile  velocity. 

The  decrease  in  the  projectile  velocity  as  it  moves  through  the  rocket  motor  is  shown  in 
Figure  5.  In  those  cases  where  far  case  layer  perforation  does  not  occur,  it  appears  that  the 
largest  velocity  drop  is  due  to  the  energy  expended  in  perforating  the  near  case  layer.  However, 
in  those  cases  where  far  case  layer  perforation  does  occur,  the  major  contributor  to  the  decrease 
in  the  projectile  velocity  is  the  fragmentation  and  motion  of  the  near  simulant  material. 

Apparently,  for  relatively  low  velocity  impacts,  most  of  the  impact  energy  is  absorbed  by  the 
perforation  of  the  near  case  layer.  However,  at  high  impact  velocities  the  creation  and 
propagation  of  the  near  propellant  simulant  material  absorbs  a  significantly  larger  portion  of  the 
initial  kinetic  energy  than  does  the  perforation  of  the  near  case  layer.  Hence,  at  impact  speeds 
above  approximately  1000  m/sec  the  largest  drop  in  projectile  velocity  occurs  after  debris  cloud 
formation.  It  is  interesting  to  note  that  the  fragmentation  of  the  far  propellant  simulant  material 
does  not  result  in  a  similar  velocity  drop.  This  fact  indicates  that  it  is  actually  the  motion  of  the 
debris  cloud,  rather  than  its  formation,  that  is  absorbing  the  majority  of  the  kinetic  energy  and 
causing  the  drop  in  the  projectile  velocity. 

Figure  6  shows  the  partitioning  of  the  initial  kinetic  energy  among  the  various  response 
modes  considered  in  this  study.  Of  immediate  note  is  the  fact  that  the  energies  associated  with  the 
delamination  of  the  near  and  far  case  layers  and  the  fragmentation  of  the  near  and  far  propellant 
simulant  material  layers  are  not  shown.  The  simple  explanation  for  this  is  that  these  response 
modes  typically  absorbed  less  than  0.5%  of  the  initial  kinetic  energy  of  the  impacting  projectile. 
Hence,  their  effects  on  the  motion  of  the  projectile  were  minimal  compared  with  the  modes 
shown.  As  expected,  Figure  6  shows  that  as  the  initial  impact  velocity  of  the  projectile  is 


increased  beyond  the  point  where  both  near  and  far  material  layers  are  perforated,  the  bulk  of  the 
kinetic  energy  remains  with  the  projectile  as  it  exits  the  far  wall  of  the  rocket  motor. 

One  area  of  concern  is  the  fact  that  the  model  is  predicting  (and  Figure  6  is  showing)  that 
as  impact  velocity  increases,  the  amount  of  energy  associated  with  the  movement  of  the  near 
debris  cloud  remains  rather  large  compared  to  the  initial  impact  energy.  It  would  be  expected  that 
as  velocity  increases,  the  projectile  would  eventually  just  perforate  the  near  propellant  simulant 
material  without  creating  a  debris  cloud  (i.e.  a  'clean'  perforation).  This  then  is  at  least  one  area  of 
the  model  that  needs  to  be  addressed  in  subsequent  studies  of  impact  energy  partitioning. 
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7.0  CONCLUSIONS  AND  RECOMMENDATIONS 

In  this  study,  a  first-principles-based  mathematical  model  was  developed  to  determine  the 
partitioning  of  the  kinetic  energy  of  an  impacting  projectile  among  various  solid  rocket  motor 
failure  modes.  Failure  modes  considered  in  the  analyses  included  near  and  far  case  perforation, 
near  and  far  case  delamination,  and  fragmentation  of  the  propellant  simulant  material.  Energies 
involved  in  material  fragmentation  were  calculated  using  a  fragmentation  scheme  based  on  a 
procedure  developed  in  a  previous  impact  study  utilizing  propellant  simulant  materials.  Two 
impact  scenarios  were  considered.  In  the  first,  the  impact  velocity  is  sufficiently  high  so  that  the 
projectile  perforates  the  near  case  and  propellant  simulant  layers,  but  not  high  enough  so  that  it 
can  perforate  the  far  case  and  simulant  material  layers.  In  the  second,  the  impact  velocity  is  high 
enough  so  that  the  projectile  is  able  to  perforate  both  sides  of  the  rocket  motor  and  exit  with 
some  residual  kinetic  energy.  The  complexity  of  the  problem  required  certain  assumptions  and 
simplifications  to  bring  the  problem  down  to  solvable  proportions. 

The  resulting  model  was  found  to  be  capable  of  predicting  a  variety  of  response 
characteristics  of  analog  solid  rocket  motors  under  high  speed  projectile  impact.  Worthy  of  note 
is  the  fact  that  the  model  correctly  predicted  the  ballistic  limit  of  the  graphite/epoxy  case  and  the 
propellant  simulant  materials.  For  the  most  part,  the  numerical  results  generated  with  the  model 
over  a  wide  range  of  impact  velocities  showed  trends  that  would  be  expected.  Naturally,  since 
the  model  was  developed  without  any  experimental  validation,  the  next  step  would  be  to  perform 
a  series  of  tests  that  would  ultimately  either  validate  the  model  or  highlight  its  weak  points. 
Experimental  test  results  are  also  required  for  the  selection  of  appropriate  values  of  the  various 
user-controlled  parameters  employed  in  the  development  of  this  model. 
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non  nnonnnnnnnonnnnnnnnnnn 


PROGRAM  I MSDSL 

IMPLICIT  DOUBLS  PRECISION  (A-H.O-Z) 

DOUBLE  PRECISION  LP , NL ,  KID , RPR , KPS , MP , MDC 
CHARACTER* 1  OPT 

OPEN ( 1 , FILE* ' IMDATA ' , STATUS- ’ UNKNOWN ' ) 
OPEN ( 2 , PILE- 1 BDELOUT ’ , STATUS- ' UNKNOWN ' ) 


READ  MATERIAL  AND  IMPACT  PARAMTERS  . . . 

RHOPR  .  PROJECTILE  DENSITY  (KG/M* 3) 

RP  .  PROJECTILE  RADIUS  (CM) 

LP  .  PROJECTILE  LENGTH  (CM) 

KPR  .  PROJECTILE  BULK  MODULUS  (K/M‘2) 

SPR  .  SLOPE  OF  PROJECTILE  US-UP  CURVE 

VP  .  IMPACT  VELOCITY  (FT/ SEC) 

NL . NUMBER  OF  LAYERS  IN  GR/EP  CASE 

TL .  LAYER  THICKNESS  IN  GR/EP  CASE 

SIGUGE  .  GR/EP  ULT  STRENGTH  (N/M‘2) 

RHOGB  .  GR/EP  DENSITY  (KG/M‘3) 

GC .  PRP  SIM  MATL  CRIT  ENERGY  REL  (N/M‘2) 

RHOPS  .  PRP  SIM  MATL  DENSITY  (KG/M‘3) 

KID .  PRP  SIM  MATL  MODE- 1  SIF  (PA-/M) 

KPS .  PRP  SIM  MATL  BULK  MODULUS  (N/M‘2) 

SPS .  PRP  SIM  MATL  US-UP  CURVE  SLOPE 

THC,THC1  ...  PRP  SIM  MATL  CONE  ANGLES  (DEG) 

AMIN .  PRP  SIM  MATL  MIN  FRAG  SIZE  (MM) 

H .  PRP  SIM  MATL  THICKNESS  (CM) 

SIGUPS  .  PRP  SIM  MATL  ULT  STRENGTH  (N/M‘2) 


READ (1,5)  RHOPR, RP,LP, KPR, SPR, VP1 

5  FORMAT ( 3F10 . S , E10 . 3 , 2F10 . 5 ) 

READ (1,6)  NL,TL, SIGUGE, RHOGE,GC 

6  FORMAT ( 2F10 . 5 , E10 . 3 , 2F10 . 5 ) 

READ (1,7)  RHOPS, KID, KPS, SPS, THC,THC1, AMIN, H, SIGUPS 

7  FORMAT(F10. 5, 2E10. 3 , F10. 5, 2F5. 2 , 2F10. 5,E10. 3 ) 

READ (1,8)  OPT 

8  FORMAT (Al) 

CONVERT  ALL  THICKNESS  TO  M,  SPEEDS  TO  M/SEC,  AND  MASSES  TO  KG 


C 

C 


C 


C 


RP-RP/ 100.0 
LP-LP/ 100.0 
VP-VP1*0 . 3048 
TL-TL/ 1000.0 
AMIN-AMIN/ 1000.0 
H-H/ 100.0 

PI-3.141592 

AH-PI*RP*RP 

MP-P I *RHOPR* RP *RP  *LP 

EINIT«MP*VP*VP/2 . 0 

TGE«NL*TL 
D-MP/ (RHOGE*AH) 

T1*D/(TGE+D) 

T2- ( 2 . 0*SIGUGE ) / ( 3 . 0*RHOGE ) 

VBL-DSQRT ( T2 * (1.0/ (T1*T1*T1)-1.0) ) 

IF  (VBL.GE.VP)  THEN 
WRITE (2, 50)  VP , VBL 

50  FORMAT ( ' THE  SPECIFIED  IMPACT  VELOCITY  (’,F6.1,’  M/S)  IS  LESS  THAN 
$THE  CALCULATED’ ,/, 'BALLISTIC  LIMIT  VELOCITY  (’,F6.1’  M/S)  OF  THE  G 
$R/EP  CASE. ’,///, ’****•  PROGRAM  TERMINATED  •****’) 

STOP 
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c 

c 


c 


c 


c 


c 

c 


END  If 

• 

EPBRFGB-MP*VBL*VBL/2 . 0 
AN-9. 0* AH 

T3-1 . 0-DEXP ( -VP/VBL ) 

EDELGE-AM*GC* ( NL- 1 . 0 ) *T3 

CPS-DSQRT ( KPS/RHOPS ) 

Dl-HP/ (RHOPS*AH) 

T4-D1/ (H+Dl) 

TS- ( 2 . 0*SIGUPS ) / ( 3 . 0*RHOPS ) 

UBL-DSQRT (T5  * (1.0/ (T4*T4*T4 ) -1 . 0 ) ) 

VF-DSQRT ( ( VP*VP+T2 ) * (T1*T1*T1 ) -T2 ) 

H1-H*DBXP ( ( 1 . 0-VF/UBL ) /4 . 0 ) 

RC- ( H-Hl ) *DTAN ( PI *THC/ 180.0) 

G1-2.0*(AMIN*AMIN*AMXN) /27 . 0 
G2-(3.0/8.0)*( (RC*RC)*(H-Hl)/3.0+(RP*RP)*Hl) 

GG— G1-G2 
HH— AMIN*  AMIN/ 9.0 
DISC«GG*GG+4 . 0*HH*HH*HH 
T6-(-GG+DSQRT(DISC) )/2.0 
T7«(-GG-DSQRT(DISC) )/2.0 
AX1-T6**{1. 0/3.0) 

AX2-T7**(1. 0/3.0) 

AMAX-AX1+AX2-AMIN/3 . 0 

ACOSF— ( 3 . 0/8 . 0 ) * (RC*RC* ( H-Hl ) /3 . 0+RP*RP*Hl ) / ( AMAX-AMIN) 

PSVOL- ( 8 . 0/3 . 0 ) *PI*ACOEF* ( AMAX-AMIN) 

T8- ( KID  *  KID ) / ( 2 . 0  *RHOPS  *CPS  *CPS ) 

T9- (8 . 0/3 . 0 ) *G2/ (AMAX-AMIN) 

T10-DSQRT ( RC*RC+ ( H-Hl ) * ( H-Hl ) ) 

EFRAGCR-T8* ( 5 . 0*PI*T9*DLOG ( AMAX/AMIN ) +PI*RC*T10 ) 

MDC— ( PI *RHOPS/3 . 0 ) *RC*RC* ( H-Hl ) +RHOPS*RP*RP*Hl 
CPR-DSQRT ( KPR/RHOPR ) 

RHORAT-RHOPS/RHOPR 

AA-SPR-SPS  *RHORAT 

BB-2 . 0*SPR*VF+CPR+CPS*RHORAT 

CC— CPR*VP+SPR*VF*VF 

DPP-BB*BB-4.0*AA*CC 

UP* ( BB-DSQRT ( DPP ) ) / ( 2 . 0 * AA ) 

EDELIV-PI* (RP*RP*RHOPS* (H/2 . 0) ) * ( 2 . 0*UP) * ( 2 . 0*UP) /2 . 0 

EFRAGMV-EDELIV-EFRAGCR 

VDC-DSQRT ( 2 . 0*EFRAGMV/MDC ) 

EPROJR-EINIT-EPERFGE-EDELGE-EFRAGCR-EFRAGMV 
VR-DSQRT ( 2 . 0*EPROJR/MP ) 

IF  (OPT.EQ. 'L' )  GOTO  100 

H11-H*DEXP( (1.0-VR/UBL) /4.0) 

RC1- ( H-Hll ) *DTAN ( PI*THC1/ 180 . 0 ) 

Gll-Gl 

G21-(3.0/8.0)*( (RCl*RC)*(H-Hll)/3.0+(RP*RP)*Hll) 

GG1— G11-G21 
HH1-HH 

DISC2»GGl*GGl+4 . 0*HH1*HH1*HH1 
T61-(-GGl+SQRT(DISC2) )/2.0 
T71-(-GGl-SQRT(DISC2) ) / 2 . 0 
AX11-T61**(1. 0/3.0) 

AX21»T71**(1. 0/3.0) 

AMAX1-AX11+AX2 l-AMIN/3 . 0 

ACOEF1«(3.0/8.0)*(RC1*RC1*(H-H11)/3.0+RP*RP*H11)/(AMAX1-AMIN) 
PSVOL1- (8 .0/3 .0) *PI*ACOEFl* ( AMAX1-AMIN) 
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c 


c 


c 


c 


c 


c 


c 


TB1-TS 

T91»{8. 0/3.0) *021/ (AMAX1-AMIN) 

T101-DSQRT (RCl*RCl+( H-Hl 1 ) * ( H-Hl 1 ) ) 

EFRGCR-T81* ( S . 0*PI*T91*DLOG (AMAX1/AMIN) +PI*RC1*T101 ) 

VI-DSQRT(VR*VR-BFRGCR/(MP/2.0) ) 

EDLGK-AM*GC* ( NL-1 . 0 ) * ( 1 . O-OBXP ( -VI /VBL ) ) 

IP  (VI.LT.VBL)  THEN 

EPROJF-O.O 

END  IP 

IP  (VI .GT. VBL)  THEN 
EPROJF-EPROJR-EFRGCR-EPERFCE-EDLGE 
END  IP 

VRF-DSQRT(EPROJF/(MP/2.0) ) 

100  WRITE (2, 101) 

101  FORMAT ( / ,  ’ IMPACT  ENERGY  PARTITIONING  FOR  CYLINDRICAL  FRAGMENT', 

S/. 'IMPACT  ON  A  FILAMENT-HOUND  CASE  WITH  PROPELLANT  SIMULANT  ...') 

RHOPR-RHOPR/ 1000 . 0 
RP-RP* 100.0 
LP»LP* 100.0 
MP-MP* 1000.0 

WRITS ( 2 , 102 )  RHOPR , KPR , SPR , CPR , RP , LP , MP , VP , VP1 

102  FORMAT ( /  ,  ' PROJECTILE  PROPERTIES  5X, ' DENSITY  ...',P10.3, 

$'  GM/CM‘3' ,/,5X, 'BULK  MODULUS  ...*,E10.3,'  N/M‘2 ',/, 5X, ' US-UP  SLOP 

SB  . . . . ' , F10.3, /,5X, 'AD  SOUND  SP  - *,F10.3,'  M/SEC' ,/, 5X, 'RADIUS 

$ . '  ,P10.3,'  CM' ,/,5X, 'LENGTH . ',P10.3,'  CM',/,SX, 

$  '  MASS  . '  ,F10.3,  '  GM’ ,/,5X, 'IMP  VEL . ',F10.3, 

S'  M/SBC  -  ' , F10.3, '  FT/SEC') 

RHOGB-RHOGE/ 1000.0 
TL«TL*1000.0 
TGE-TGE* 1000.0 

WRITE (2,103)  RHOGE , NL , TL, TGE , SIGUGE , GC , VBL 

103  FORMAT ( 'GR/EP  CASE  PROPERTIES  ...',/, 5X, 'DENSITY  . ', 

SF10.3,'  GM/CM‘3' ,/,5X, '#  LAYERS  . '  ,  F10. 3 ,  /  ,  5X, 'LAYER  T 

SHKNBSS . '  ,  F10.3,  '  MM' ,/,  5X, 'TOTAL  THKNESS  . ',F10.3, 

$'  MM' ,/,5X, 'ULT  STRENGTH  . ',E10.3,'  N/M‘2 ’,/, 5X, 'CRIT  ENER 

$  REL  RATE  ...',F10.3,'  N/M‘2 ',/, 5X, ' BALL  LIM  VEL  . ',F10.3, 

$'  M/SEC') 

RHOPS-RHOPS/ 1000 . 0 
AMIN-AMIN* 1000.0 
H«H*100.0 

WRITE (2,104)  RHOPS , KID ,KPS,SPS,CPS,S IGUPS , AMIN , H , UBL 

104  FORMAT ( 'PROPELLANT  SIMULANT  MATERIAL  PROPERTIES  ...',/, 5X, 

$ 'DENSITY  . ' ,F10.3, '  GM/CM‘3 ’,/, 5X, 'DYN  FRACT  STR  ...', 

$E10.3, '  N/M*2' ,/,5X, 'BULK  MODULUS  - ',E10.3,'  N/M‘2 ’,/, 5X, ' US-UP 

$  SLOPE  . ' ,F10.3,/,SX, 'AD  SOUND  SP  . ',F10.3,'  M/SEC', /,5X, 

$'ULT  STRENGTH  _ ',610.3,'  N/M‘2 ’,/, 5X, 'MIN  FRAG  SIZE  ...', 

$F10. 3, ’  MM' ,/,5X,  "THICKNESS  . ',F10.3,'  CM* , / , 5X, ’ BALL  LIM  VE 

$L - ’  , F10.3,  '  M/SEC’  ) 

H1«H1*100.0 

HC»H-H1 

RC-RC* 100.0 

MDC»MDC* 1000.0 

AMAX>AMAX* 100.0 

ACOEF*ACOEF* 100. 0*100.0 

PSVOL«PSVOL* 100 . 0* 100 .0*100.0 

WRITE (2,105)  EINIT, EPERFGE , VF , EDELGE , EDELIV , EFRAGCR , ACOEF , AMAX , 

$RC , HC , THC , PS VOL , MDC , UP , EFRAGMV , VDC , EPROJR , VR 
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c 


c 


c 


105  FORMAT ( / , ' ENERGY  PARTITIONING  DUS  TO  FIRST  HALL  PERFORATION  ...',/ 

$,  5X,  '  INITIAL  PROJECTILE  KINETIC  ENERGY  . '.F10.3,'  J',/ 

$,  5X,  *  ENERGY  REQD  TO  PERF  GR/EP  CASE  . ',F10.3,'  J*,/, 

$10X, 'RES  PRJ  VEL  GR/BP  PERF . *,F10.3,'  M/SEC',/ 

$,  5X,  'ENERGY  ABS  IN  GR/EP  DELAMINATION  . ',F10.3,'  J',/ 

$ ,  5X ,  '  ENERGY  DELIVERED  TO  PRP  SIM  MATL . ',F10.3,'  J',/, 

$ ,  5X ,  '  ENERGY  REQD  TO  FRAG  PRP  SIM  MATL  . ',F10.3,'  J',/, 

$10X, 'CUM  FRAG  FCN  COEFFICIENT  . ',F10.3,'  CM‘2',/, 

$10X,  'MAX  PRP  SIM  MATL  FRAG  DIA - *,F10.3,'  CM',/, 

$10X, ' PRP  SIM  MATL  CONE  PLUG  RAD  ...',F10.3,'  CM',/, 

S10X, 'PRP  SIM  MATL  CONE  PLUG  HT  _ ',F10.3,'  CM’,/, 

$10X, 'PRP  SIM  MATL  CONE  1/2-ANG  - *,F10.3,'  DEG',/, 

$10X, ' PRP  SIM  MATL  FRAG  VOLUME  . ',F10.3,'  CM‘3',/, 

$10X,'PRP  SIM  MATL  DEB  CLD  MASS  _ *,F10.3,'  GM',/, 

$10X, 'PRP  SIM  MATL  PARTICLE  VEL  _ ',F10.3,'  M/SEC',/, 

5 5 X, 'ENERGY  AVAIL  TO  MV  PRP  SIM  MATL  DEB  CLD  - ',F10.3,'  J',/, 

$10X, 'PRP  SIM  MATL  DEB  CLD  VEL  . *,FI0.3,'  M/SEC',/ 

$,5X, 'RESIDUAL  PROJECTILE  KINETIC  ENERGY  . ',F10.3,'  J’,/, 

$10x,  ' RESID  PROJ  VELOCITY . *,FI0.3,'  M/SEC') 

Fl-EPERFGE/EINIT 

F2-EDELGE/EINIT 

F3-EFRAGCR/EINIT 

F4«EFRAGMV/EINIT 

F5-EPROJR/EINIT 

WRITE(2, 106)  FI, F2,F3,F4,F5 

106  FORMAT (/, ' EPERFGE/EINIT  »  ' , F10 . 3 , / , * EDELGE/EINIT  *  ',F10.3, 

$/, 'EFRAGCR/EINIT  *  ' , F10 . 3 , / , ' EFRAGMV/EINIT  *  ',F10.3,/, 

$ ' EPROJR/EINIT  -  * , F10 . 3 ) 

IF  (OPT.EQ. 'L')  THEN 
FS-F1+F2+F3+F4+FS 
WRITE (2 ,107)  FS 

107  FORMAT ( ' SUM  OF  RATIOS  *  ',F10.3,/) 

IF  (VR.GT.VBL)  WRITE(2,108) 

108  FORMAT ('*****  WARNING  *****  RESID  PROJ  VEL  EXCEEDS  BALLISTIC  LIMIT 
$  OF',/, 'REAR  GR/EP  CASE  WALL.  RERUN  USING  HIGH  VELOCITY  OPTION.') 

GOTO  200 
END  IF 

H11=*H11*  100.0 

HC1-H-H11 

RC1*RC1* 100.0 

AMAX1=AMAX1* 100.0 

ACOEF1«ACOEF1*100 . 0*100 . 0 

PSVOLl^PSVOLl* 100 .0*100.0*100.0 

WRITE (2,150)  EPROJR, EFRGCR , ACOEF1 , AMAX1 , RC1 , HC1 , THC1 , PSVOL1 , VI , 
$EDLGE , EPERFGE , EPROJF , VRF 

150  FORMAT (/, 'ENERGY  PARTITIONING  DUE  TO  SECOND  WALL  PERFORATION  ...', 

$/,5X, 'RESID  PROJ  KIN  ENERGY  AFTER  1ST  WALL  PERF  _ ',F10.3,'  J', 

$/,5X, 'ENERGY  REQD  TO  FRAG  PRP  SIM  MATL . ’,F10.3,’  J',/, 

$10X, ' PRP  SIM  MATL  FRAG  FCN  COEFF  . ’,F10.3,'  CM‘2',/, 

$10X, 'MAX  PRP  SIM  MATL  FRAG  DIA  . ',F10.3,'  J',/, 

$10X, ' PRP  SIM  MATL  PLUG  CONE  RAD . ',F10.3,'  CM',/, 

$10X,  '  PRP  SIM  MATL  PLUG  CONE  HT . ',F10.3,'  CM',/ 

$10X, ' PRP  SIM  MATL  CONE  1/2-ANG  . ',F10.3,'  DEG',/, 

$10X,  ' PRP  SIM  MATL  FRAG  VOLUME  . ',F10.3,'  CM‘3',/, 

$10X, 'POST  PRP  SIM  MATL  FRAG  PROJ  VEL  ...’,F10.3,’  M/SEC',/ 

$,5X, 'ENERGY  ABS  BY  IN  GR/EP  DELAMINATION  . ’,F10.3,'  J',/ 

$,5X, 'ENERGY  REQD  TO  PERF  GR/EP  CASE  . ',F10.3,'  J',/ 

$,5X, 'FINAL  RESIDUAL  PROJECTILE  KINETIC  ENERGY  . ',F10.3,'  J',/, 

$10X, 'FINAL  RES  PROJ  VEL . ’,F10.3,'  M/SEC') 

F6-EFRGCR/EINIT 
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P7-KDLGI/XXNXT 
F8-BPERFGB/EINIT 
F9-SPR0JF/EINIT 
FS-F1+F2+F3+F4+F6+F7+F8+F9 
WRITE(2, 151)  F6,F7,F8,F9,FS 

151  FORMAT(/, 'EFROCR/BINIT  -  '  ,  F10 . 3 ,  / ,  '  EOLGE/EINIT  *  *,F10.3,/, 

$ ' EPERFGB/EINXT  »  1 ,F10.3, /, ’EPROJF/EINIT  -  ' , F10. 3, / , • SUM  OF  RATI 
$OS  -  ’ , F10. 3 ) 

C 

200  STOP 
END 
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APPENDIX  B  -  SAMPLE  INPUT  AND  OUTPUT  FILES 
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Sample  Input  File  IMDATA 


7830.0 

0.5 

2.0 

1.635E+11 

1.55 

11000.0 

6.0 

0.4 

0.200E+10 

1530.0 

400.0 

1780.0 

0.128E+07 

0.817E+10 

1.95 

45.0  45.0 

0.1 

2.54  0.590E+07 
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Sample  Output  File  for  a  Low  Velocity  Impact ... 


IMPACT  ENERGY  PARTITIONING  FOR  CYLINDRICAL  FRAGMENT 
IMPACT  ON  A  FILAMENT-MOUND  CASE  WITH  PROPELLANT  SIMULANT  . . . 


rnuu x ls  rnursiu xao  ... 

DENSITY  ...  7.830  GM/CM‘3 

BULK  MODULUS  ...  0.164E+12  N/M*2 

US-UP  SLOPE  -  1.550 

AD  SOUND  SP  -  4569.598  M/SEC 

RADIUS  .  0.500  CM 

LENGTH  .  2.000  CM 

MASS .  12.299  GM 

IMP  VSL  .  304.800  M/SEC  *  1000.000  FT/SEC 

GR/EP  CASE  PROPERTIES  ... 

DENSITY  .  1.530  GM/CM‘3 

/  LAYERS .  6.000 

LAYER  THKNBSS .  0.400  MM 

TOTAL  THKNBSS .  2.400  MM 

ULT  STRENGTH  .  0.200E+10  N/M*2 

CRIT  ENER  REL  RATE  ...  400.000  N/M*2 

BALL  LIM  VEL .  250.502  M/SEC 


PROPELLANT  SIMULANT 

DENSITY  . 

DYN  FRACT  STR 
BULK  MODULUS  . 
US-UP  SLOPE  . . 
AD  SOUND  SP  .. 
ULT  STRENGTH  . 
MIN  FRAG  SIZE 
THICKNESS  .... 
BALL  LIM  VEL  . 


MATERIAL  PROPERTIES  ... 
..  1.780  GM/CM‘3 

..  0.128E+07  N/M‘2 
..  0.817E+10  N/M‘2 
1.950 

..  2142.402  M/SEC 

..  0.590E+07  N/M‘2 
..  0.100  MM 

2.540  CM 

..  50.196  M/SEC 


ENERGY  PARTITIONING  DUE  TO  FIRST  HALL  PERFORATION  . . . 

INITIAL  PROJECTILE  KINETIC  ENERGY  .  571.323 

ENERGY  REQD  TO  PERF  GR/EP  CASE  .  385.900 

RES  PRJ  VEL  GR/EP  PERF  .  167.709  M/SEC 

ENERGY  ABS  IN  GR/EP  DELAMINATION  .  0.995 

ENERGY  DELIVERED  TO  PRP  SIM  MATL .  79.568 

ENERGY  REQD  TO  FRAG  PRP  SIM  MATL .  0.884 

CUM  FRAG  FCN  COEFFICIENT  .  0.468  CM* 2 

MAX  PRP  SIM  MATL  FRAG  DIA  -  0.674  CM 

PRP  SIM  MATL  CONE  PLUG  RAD  ...  1.125  CM 

PRP  SIM  MATL  CONE  PLUG  HT  -  1.125  CM 

PRP  SIM  MATL  CONE  1/2-ANG  _  45.000  DEG 

PRP  SIM  MATL  FRAG  VOLUME  .  2.603  CM *3 

PRP  SIM  MATL  DEB  CLD  MASS  _  3.286  GM 

PRP  SIM  MATL  PARTICLE  VEL  -  149.691  M/SEC 

ENERGY  AVAIL  TO  MV  PRP  SIM  MATL  DEB  CLD  -  78.684 

PRP  SIM  MATL  DEB  CLD  VEL  .  218.842  M/SEC 

RESIDUAL  PROJECTILE  KINETIC  ENERGY  .  104.860 

RESID  PROJ  VELOCITY .  130.581  M/SEC 


J 

J 

J 

J 

J 


J 

J 


EPERFGE/EINIT  ■ 

0.675 

EDELGE/EINIT  * 

0.002 

EFRAGCR/EINIT  - 

0.002 

EFRAGMV /EINIT  - 

0.138 

EPROJR/BINIT  - 

0.184 

SUM  OF  RATIOS  - 

1.000 

Sample  Output  File  for  a  Medium  Velocity  Impact ... 

IMPACT  ENERGY  PARTITIONING  FOR  CYLINDRICAL  FRAGMENT 
IMPACT  ON  A  FILAMENT-HOUND  CASE  WITH  PROPELLANT  SIMULANT  . . . 


PROJECTILE  PROPERTIES 


DENSITY  ... 
BULK  MODULUS 
US-UP  SLOPE  . 
AD  SOUND  SP  . 

RADIUS  . 

LENGTH  . 

MASS  . 

IMP  VEL  . 


7.830  GM/CM‘3 
..  0 . 164E+12  N/M*2 
1.5S0 

..  4569. 598  M/SEC 

. .  0.500  CM 

..  2.000  CM 

..  12.299  GM 

..  609.600  M/SEC 


2000.000  FT/SEC 


GR/BP  CASE  PROPERTIES 

DENSITY  . 

#  LAYERS  . 

LAYER  THKNESS  ... 
TOTAL  THKNESS  ... 


1.530  GM/CM*3 
6.000 
0.400  MM 
2.400  MM 


ULT  STRENGTH  .  0.200E+10  N/M*2 

CRIT  ENER  REL  RATE  ...  400.000  N/M‘2 

BALL  LIM  VEL .  250.502  M/SEC 

PROPELLANT  SIMULANT  MATERIAL  PROPERTIES  ... 


DENSITY  . 

DYN  FRACT  STR 
BULK  MODULUS  . 
US-UP  SLOPE  .. 
AD  SOUND  SP  .. 
ULT  STRENGTH  . 
MIN  FRAG  SIZE 
THICKNESS  .... 
BALL  LIM  VEL  . 


1.780 
0.128E+07 
0.817E+10 
1.950 
2142.402 
0. 590E+07 
0.100 
2.540 
50.196 


GM/CM‘3 

N/M‘2 

N/M*2 

M/SEC 

N/M‘2 

MM 

CM 

M/SEC 


ENERGY  PARTITIONING  DUE  TO  FIRST  WALL  PERFORATION 

INITIAL  PROJECTILE  KINETIC  ENERGY  . 

ENERGY  REQD  TO  PERF  GR/EP  CASE  .... 

RES  PRJ  VEL  GR/EP  PERF  . 

ENERGY  ABS  IN  GR/EP  DELAMINATION  . 

ENERGY  DELIVERED  TO  PRP  SIM  MATL  . 

ENERGY  REQD  TO  FRAG  PRP  SIM  MATL  . 


J 

J 


536.763 


2285.291 
385.900 
M/SEC 

1.290  J 
775.711  J 
2.969  J 


CUM 

FRAG  FCN 

COEFFICIENT  .. 

•  •  • 

1.367 

CM*  2 

MAX 

PRP 

SIM  MATL  FRAG  DIA  . 

•  •  • 

1.159 

CM 

PRP 

SIM 

MATL 

CONE  PLUG  RAD 

•  •  • 

2.315 

CM 

PRP 

SIM 

MATL 

CONE  PLUG  HT  . 

•  •  • 

2.315 

CM 

PRP 

SIM 

MATL 

CONE  1/2-ANG  . 

•  •  • 

45.000 

DEG 

PRP 

SIM 

MATL 

FRAG  VOLUME  . . 

•  •  • 

13.167 

CM' 3 

PRP 

SIM 

MATL 

DEB  CLD  MASS  . 

•  •  • 

23.223 

GM 

PRP 

SIM 

MATL 

PARTICLE  VEL  . 

•  •  • 

467.389 

M/SEC 

ENERGY  AVAIL 

TO  MV  PRP  SIM  MATL 

DEB 

OLD  • • • • 

772. 

PRP 

SIM 

MATL 

DEB  CLD  VEL  . . 

•  •  • 

257.973 

M/SEC 

RESIDUAL  PROJECTILE  KINETIC  ENERGY  .  1122.391 

RESID  PROJ  VELOCITY  .  427.215  M/SEC 


EPERFGE/EINIT  - 

0.169 

EDELGE/EINIT  = 

0.001 

EFRAGCR/EINIT  - 

0.001 

BFRAGMV/EINIT  » 

0.338 

EPROJR/EINIT  - 

0.491 

40 


ENERGY  PARTITIONING  DUB  TO  SECOND  WALL  PERFORATION  ... 

RES ID  PROJ  KIN  ENERGY  AFTER  1ST  WALL  PERF  -  1122.391  J 

ENERGY  RBQD  TO  FRAG  PRP  SIM  MATL .  2.682  J 

PRP  SIM  MATL  FRAG  FCN  COEFF .  1.166  CM‘2 

MAX  PRP  SIM  MATL  FRAG  DIA .  1.109  J 

PRP  SIM  MATL  PLUG  CONE  RAD .  2.152  CM 

PRP  SIM  MATL  PLUG  CONE  HT .  2.152  CM 

PRP  SIM  MATL  CONE  1/2-ANG .  45.000  DEG 

PRP  SIM  MATL  FRAG  VOLUME  .  10.735  CM*  3 

POST  PRP  SIM  MATL  FRAG  PROJ  VEL  ...  426.704  M/SEC 

ENERGY  ABS  BY  IN  GR/BP  DELAMINATION  .  1.156  J 

ENERGY  REQD  TO  PERF  GR/EP  CASE  .  385.900  J 

FINAL  RESIDUAL  PROJECTILE  KINETIC  ENERGY  .  732.652  J 

FINAL  RES  PROJ  VEL .  345.162  M/SEC 

EFRGCR/EINIT  -  0.001 

EDLGB/EINIT  -  0.001 

EPBRFGE/BINIT  -  0.169 

EPROJF/EINIT  -  0.321 

SUM  OF  RATIOS  -  1.000 


r 
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Sample  Output  File  for  a  High  Vdocity  Impact ... 

IMPACT  ENERGY  PARTITIONING  FOR  CYLINDRICAL  FRAGMENT 
IMPACT  ON  A  FILAMENT-WOUND  CASE  WITH  PROPELLANT  SIMULANT  ... 


rnwiivuiis  rnurinuu  ... 

DENSITY  ...  7.830  GM/CM‘3 

BULK  MODULUS  ...  0.164B+12  N/M‘2 

US-UP  SLOPE  -  1.550 

AD  SOUND  SP  -  4569. S98  M/SEC 

RADIUS  .  0.500  CM 

LENGTH  .  2.000  CM 

MASS  .  12.299  GM 

IMP  VEL  .  3352.800  M/SEC  •  11000.000  FT/SEC 

GR/EP  CASE  PROPERTIES  ... 

DENSITY  .  1.530  GM/CM‘3 

#  LAYERS .  6.000 

LAYER  THKNESS .  0.400  MM 

TOTAL  THKNESS  .  2.400  MM 

ULT  STRENGTH  .  0.200E+10  N/M‘2 

CRIT  ENER  REL  RATE  ...  400.000  N/M‘2 

BALL  LIM  VEL .  250.502  M/SEC 


PROPELLANT  SIMULANT  MATERIAL  PROPERTIES 


DENSITY  . 

DYN  FRACT  STR 
BULK  MODULUS  . 
US-UP  SLOPE  .. 
AD  SOUND  SP  .. 
ULT  STRENGTH  . 
MIN  FRAG  SIZE 
THICKNESS  .... 
BALL  LIM  VEL  . 


1.780 
0. 128E+07 
0.817E+10 
1.950 
2142.402 
0. 590E+07 
0.100 
2.540 
50.196 


GM/CM' 

N/M‘2 

N/M‘2 

M/SEC 

N/M‘2 

MM 

CM 

M/SEC 


ENERGY  PARTITIONING  DUE  TO  FIRST  WALL  PERFORATION  . . . 

INITIAL  PROJECTILE  KINETIC  ENERGY  .  69130.046 

ENERGY  REQD  TO  PERF  GR/EP  CASE  .  385.900 

RES  PRJ  VEL  GR/EP  PERF  .  3229.187  M/SEC 

ENERGY  ABS  IN  GR/EP  DELAMINATION  .  1.414 

ENERGY  DELIVERED  TO  PRP  SIM  MATL .  22508.139 

ENERGY  REQD  TO  FRAG  PRP  SIM  MATL .  3.602 

CUM  FRAG  FCN  COEFFICIENT  .  1.630  CM*2 

MAX  PRP  SIM  MATL  FRAG  DIA  ....  1.267  CM 

PRP  SIM  MATL  CONE  PLUG  RAD  ...  2 . 540  CM 

PRP  SIM  MATL  CONE  PLUG  HT  ....  2.540  CM 

PRP  SIM  MATL  CONE  1/2-ANG  ....  45.000  DEG 

PRP  SIM  MATL  FRAG  VOLUME  .  17.160  CM* 3 

PRP  SIM  MATL  DEB  CLD  MASS  _  30.546  GM 

PRP  SIM  MATL  PARTICLE  VEL  _  2517.666  M/SEC 

ENERGY  AVAIL  TO  MV  PRP  SIM  MATL  DEB  CLD  -  22504.537 

PRP  SIM  MATL  DEB  CLD  VEL  .  1213.879  M/SEC 

RESIDUAL  PROJECTILE  KINETIC  ENERGY  .  46234.593 

RESID  PROJ  VELOCITY  .  2741.939  M/SEC 


J 

J 

J 

J 

J 


J 

J 


EPERFGE/EINIT  »  0.006 
BDELGE/EINIT  *  0.000 
EFRAGCR/EINIT  =  0.000 
EFRAGMV/EINIT  *  0.326 
EPROJR/EINIT  ■  0.669 
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BN* ROY  PARTITIONING  DUB  TO  SECOND  WALL  PERFORATION  . . . 

RBSID  PROJ  KIN  BNBRGY  AFTER  1ST  WALL  PBRF  _  46234.593  J 


ENERGY  RBQD  TO  FRAG  PRP  SIM  MATL .  3.602  J 

PRP  SIM  MATL  FRAG  FCN  COEFF .  1.630  CM'2 

MAX  PRP  SIM  MATL  FRAG  DIA .  1.267  J 

PRP  SIM  MATL  PLUG  CONE  RAD .  2.540  CM 

PRP  SIM  MATL  PLUG  CONE  HT .  2.540  CM 

PRP  SIM  MATL  CONE  1/2-ANG .  45.000  DEG 

PRP  SIM  MATL  FRAG  VOLUME  .  17.160  CM~3 

POST  PRP  SIM  MATL  FRAG  PROJ  VEL  ...  2741.832  M/SEC 

ENERGY  ABS  BY  IN  GR/EP  DELAMINATION  .  1.414  J 

ENERGY  REQD  TO  PERF  GR/EP  CASE  .  385.900  J 

FINAL  RESIDUAL  PROJECTILE  KINETIC  ENERGY  .  45843.678  J 


FINAL  RES  PROJ  VEL 


2730.323  M/SEC 


EFRGCR/EINIT  *  0.000 
EDLGE/BINIT  -  0.000 
EPBRFGB/EINIT  »  0.006 
EPROJF/EINIT  «  0.663 
SUM  OF  RATIOS  «  1.000 
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Generic  Solid  Rocket  Motor  Impacted  by  a  Cylindrical  Projectile 
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Figure  2.  Flat  Plate  Model  of  Solid  Rocket  Motor 


Propellant  Simulant  Material  Cavity  Model 
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Figure  4  Propellant  Simulant  Debris  Cloud  and  Final  Projectile  Velocity  vs  Initial 
Impact  Velocity 
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Velocity  Value  Station  (1... Impact  Velocity,  2... After  Primary  Case  Perforation,  3... After 
Primary  Propellant  Simulant  Fragmentation,  4... After  Secondary  Propellant  Simulant 
Fragmentation,  5. ..After  Secondary  Case  Perforation 
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Figure  6.  Partitioning  of  Initial  Projectile  Kinetic  Energy 

(NCL  Near  Case  Layer,  FCL  ...  Far  Case  Layer) 
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